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Abstract We have developed the Monte Carlo simulation program Jewel 1.0 (Jet Evolution With Energy 
Loss), which interfaces a perturbative final state parton shower with medium effects occurring in ultra- 
relativistic heavy ion collisions. This is done by comparing for each jet fragment the probability of further 
perturbative splitting with the density-dependent probability of scattering with the medium. A simple 
hadronisation mechanism is included. In the absence of medium effects, we validate Jewel against a set 
of benchmark jet measurements. For elastic interactions with the medium, we characterise not only the 
medium-induced modification of the jet, but also the jet-induced modification of the medium. Our main 
physics result is the observation that coUisional and radiative medium modifications lead to characteristic 
differences in the jet fragmentation pattern, which persist above a soft background cut. We argue that this 
should allow to disentangle coUisional and radiative parton energy loss mechanisms by measuring the n-jet 
fraction or a class of jet shape observables. 



1 Introduction 

In ultra-relativistic heavy ion collisions, the produced QCD 
matter reduces significantly the energy of high transverse 
momentum partons. Experiments at the Relativistic Heavy 
Ion Collider RHIC give strong support to this picture [1- 
4] . Most of the experimental evidence comes from studying 
the leading hadronic fragments of the parent partons via 
single inclusive hadron spectra and jet-like particle cor- 
relations. Since the medium-induced suppression of these 
spectra persists unattenuated up to the highest transverse 
momenta accessed experimentally at RHIC, it is likely to 
play a dominant role also in the much wider transverse 
momentum range soon to be explored in heavy ion colli- 
sions at the CERN Large Hadron Colhder LHC [5-7]. 

On general grounds, the energy lost by the leading 
fragment in a parton shower must manifest itself in the as- 
sociated particle yield. So, parton energy loss is expected 
to result in a medium modification of the entire parton 
fragmentation pattern ('jet quenching') and not only in 
the suppression of leading hadrons. The corresponding 
study of medium-modified jets beyond their leading frag- 
ments is of great interest for several reasons, in particular: 
i) At the LHC a larger fraction of the entire medium- 
modified jet fragmentation pattern will become accessible 
above background, ii) Studying the distribution of sub- 
leading fragments is likely to discriminate between dif- 
ferent microscopic mechanisms conjectured to under ly jet 
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quenching, thereby helping to characterise more precisely 
the properties of matter tested by jet quenching, iii) Mod- 
elling the distribution of subleading jet fragments is essen- 
tial for any operational procedure aiming at disentangling 
jets from background or characterising the jet-induced 
modification of the background. Such reasons motivate 
the development of tools which account dynamically for 
the interaction between jet and medium, and which model 
medium-modified jets on the level of multi-particle final 
states. 

The interactions of a partonic jet component with some 
target (typically, a quark or a gluon) in the medium can be 
elastic or inelastic. In the elastic case, the projectile parton 
can lose energy only by transferring recoil momentum to 
the target ('coUisional parton energy loss' [8-16]). In the 
inelastic case, the dominant source of energy degradation 
of the partonic projectile is the medium-induced break- 
up of the projectile via gluon radiation ('radiative parton 
energy loss' [17-22]), resulting in two or more projectile 
components of lower energy. 

The Monte Carlo (MC) technique provides a powerful 
tool for the simulation of multi-particle final states. In the 
absence of medium-effects, the perturbative dynamics of 
parton fragmentation maps with known accuracy onto a 
probabilistic iteration of parton splittings, and the Monte 
Carlo technique is widely used for the simulation of final 
state parton showers in the vacuum [23-25]. In the pres- 
ence of medium effects, we expect that a parton shower 
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can contribute to understanding 'jet quenching', in par- 
ticular for the following reasons: 

1. Baseline for the study of medium-modified multi-particle 
states. 

The standard experimental procedure of characteris- 
ing medium-modifications of jet-like multi-particle fi- 
nal states is to compare to an experimental baseline in 
which medium effects are absent. A realistic modelling 
should be able to parallel this experimental procedure. 

2. Exact implementation of energy and momentum con- 
servation. 

Multi-particle final states are known to be sensitive to 
trigger biases and they are constrained to share the 
energy and momentum of a parent parton. For the de- 
scription of this constrained dynamics, exact energy- 
momentum conservation is likely to be important. It 
may be numerically more important than the treat- 
ment of interference effects, on which analytical calcu- 
lations of radiative parton energy loss tend to focus. 

3. Testing a wide range of microscopic mechanisms for 
the interaction between projectile and target. 

A Monte Carlo tool should allow the user to change 
easily between different descriptions of the medium 
and of the interaction between the medium and the 
partonic projectile. This is important to constrain the 
microscopic mechanism underlying parton energy loss 
and to better identify the specific properties of the 
medium tested by jet quenching. 

4. Interface with experiment. 

The modelling of realistic multi-hadron final states is 
clearly beneficial for comparing theory and data, in- 
cluding detailed detector acceptance, resolution and 
response. 

We note that the application of Monte Carlo tech- 
niques for the modelling of medium-effects leaves open 
conceptual issues in the current description of jet quench- 
ing. In particular, medium-induced gluon radiation is ex- 
pected to arise from an interference phenomenon whose 
mapping onto a probabilistic description of multi-gluon 
emission is not known. Also, while the iteration of coUi- 
sional mechanisms may be probabilistic, its interface with 
a standard parton shower involves modelling assumptions. 
In the present work, we have bypassed these fundamen- 
tal issues in order to arrive at a tool which realises the 
practically important points listed above. 

The working name of this Monte Carlo code is Jewel, 
standing for Jet Evolution With Energy Loss. In the ab- 
sence of a medium, the Jewel parton shower is ordered in 
virtuality, and strict angular ordering is enforced by phase 
space constraints. The medium is modeled as a collection 
of scattering centres whose cross sections and distribution 
in phase space can, in principle, be chosen freely. To de- 
cide whether a parton splits as in the vacuum or whether 
it interacts with a target parton in the medium, the code 
compares the mean free path between scattering centres 
with the mean lifetime of the virtual parton. Within this 
framework, Jewel aims at realising the objectives listed 
above. 



One approach currently used to connect parton en- 
ergy loss calculations to measurable quantities is based 
on energy loss probabilities ('quenching weights' [26,27]), 
defined by the probabilistic iteration of single medium- 
induced gluon emissions. Alternatively, the interpretation 
of medium-induced gluon radiation as a modification of 
parton splitting functions has motivated attempts to evolve 
partons via medium-modified DGLAP evolution equations. 
Both approaches lead to comparable results for leading 
hadronic fragments [28]. Also, coUisional energy loss mech- 
anisms have been used recently to define quenching weights 
and medium-modified DGLAP evolutions [16]. We note 
that these approaches are largely limited to leading hadron 
spectra, or special classes of two-particle correlation func- 
tions [29], mainly since they treat subleading fragments 
with kinematic approximations. Recently, there have been 
many efforts to go beyond leading hadron spectra within 
the framework of the models discussed above. This in- 
cluded works on jet broadening [30], the angular depen- 
dence of two-particle correlations [31], the effect of di- 
rected momentum transfer (a.k.a. flow) on the jet frag- 
mentation [32,33], the modification of jet multiplicity dis- 
tributions [30,34] and jet hadrochemistry [35]. We expect 
that Monte Carlo models of jet quenching will greatly ad- 
vance this field of study, mainly because they can imple- 
ment the four important features listed above. 

In recent years, several Monte Carlo models imple- 
menting parton energy loss have been developed already. 
In particular, the code HiJiNG [37] simulates complete 
nucleus-nucleus collisions and provides a simplified model 
for radiative energy loss derived from analytical calcula- 
tions. There are also codes, which simulate jet quench- 
ing without modeling the entire nucleus-nucleus collision, 
such as Pqm [38], which is a Monte Carlo implementa- 
tion of the BDMPS quenching weights, and Pyquen [36], 
which modifies the standard Pythia 6.2 jet events es- 
sentially by reducing the energy of partons in the shower 
and by adding additional gluons to that shower according 
to distributions motivated by parton energy loss calcula- 
tions. There are also two models of medium modified par- 
ton showers, namely Q-Pythia [39] with modified split- 
ting functions including induced gluon radiation in the 
BDMPS model and a model where the virtuality of par- 
tons in the shower is increased due to interactions with the 
medium thus stimulating additional gluon emissions [40]. 
All these approaches go beyond leading hadron spectra. 
In developing Jewel, our main focus was to arrive at a 
code which allows to study in detail the dynamics relat- 
ing the evolution of the parton shower to the microscopic 
modeling of the medium interactions. In particular, we 
are interested in specifying dynamically the physically in- 
teresting (but model-dependent) relation between parton 
energy loss, p^-broadening, recoil momentum, change in 
jet multiplicity, and other characteristic features of parton 
energy loss. 

This paper describes the physics encoded in Jewel 
version 1.0, and demonstrates its use for the study of jet 
quenching phenomena on the level of multi-hadron final 
states. The paper is organised as follows: In section [21 we 
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discuss the Monte Carlo final state parton shower and its 
interface with a hadronisation mechanism in the absence 
of medium effects. To arrive at a manageable inclusion 
of medium-effects on the hadron level, the hadronisation 
model used by Jewel is less sophisticated than what is 
implemented in state of the art Monte Carlo event gener- 
ators. We discuss these differences and show that Jewel 
provides an adequate description of many jet observables, 
including event shapes, n-jet fractions and intra-jet multi- 
plicity distributions. In section [3l we discuss how scatter- 
ing in the plasma is included in Jewel. Most of our dis- 
cussion will focus on the case of elastic scatterings. For this 
case, we test in section [3] a well-studied theoretical base- 
line, namely the in-medium propagation of a parton in the 
absence of parton splitting. Section |4] addresses the ques- 
tion to what extent the medium-modifications included 
in Jewel affect experimentally accessible jet observables 
in heavy ion collisions. Jewel, as any newly developed 
Monte Carlo event generator, is a work in progress. In the 
conclusions, we summarise what we have learned from the 
present exploratory study and highlight some directions 
for future work. 



2 The Monte Carlo Model in the absence of 
medium effects 

In this section, we introduce the baseline on top of which 
medium-effects are included and show that it basically re- 
produces the observed QCD radiation physics of jets in 
the vacuum. The evolution variable of a parton shower is 
not unique. We have decided to use the virtuality Q^. As 
explained in more detail in section [3l this has the advan- 
tage that the evolution variable traces the lifetime 1/Q 
of the virtual states, which facilitates the embedding of 
the parton shower in the spatiotemporal geometry of a 
medium. We interface this parton shower with a hadro- 
nisation scheme which implements the idea that colour 
neutralisation occurs locally during hadronisation. How- 
ever, the scheme invoked here is less sophisticated than 
the hadronisation prescriptions used in modern event gen- 
erators, in particular in that it does not require knowledge 
about the event-specific colour flow in the parton shower. 
This is a technical simplification, which - in contrast to 
standard treatments - allows for a straightforward exten- 
sion of the hadronisation mechanism in the presence of a 
medium. 



2.1 Final state parton shower in the absence of 
medium effects 

2.1.1 Parton evolution 

We want to describe the evolution of a parton of initial en- 
ergy E, produced in a hard scattering process. This parton 
fragments into a multi-parton final state. In the absence of 
a medium, the Jewel parton shower is closely related to 
the mass-ordered shower in the Pythia 6.4 event genera- 
tor [23]. The kinematics of each a b+c parton branching 



is given in terms of the virtuality of the parent parton and 
the momentum fraction z carried by one of its daughters. 
The probability that no splitting occurs between an initial 
and final virtuality Qi and Q/, respectively, is described 
by the Sudakov factor 



(1) 



exp 



Qf ^+(Q'^B) 



2tt 



b,c 



Here, Pa,^hc{z) are the standard LO parton splitting func- 
tions for quarks and gluons {a,b,c S {q,g}). In order to 
regularise the integral one has to define an infrared cut- 
off scale below which splittings are considered to be not 
resolvable. In practice, we require a minimal virtual mass 
Qo/2 for the daughters. Given this cut-off, the condition 
fc^ z(l — z)Q^ > ^QCD translates directly into the al- 
lowed z range 
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To avoid the Landau pole in equation ([1]) , we increase the 
nominal cut on k± by 10%. These are the same prescrip- 
tions as used in the mass-ordered Pythia cascade. With 
the no-splitting probability equation ([T]) the probability 
density Sii{Qf, Q^) for a splitting to happen at virtuality 
is given by 



^a(gf,Q') 



where 



VFa-.bc(Q' 



d(lnQ2) 



= S,iQt,Q') 



b,c 



(3) 
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dz 



a,(z(l-z)Q2) 
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Pa-^bc(^) (4) 



is the differential probability for the splitting a ^ & + c at 
Q^, and the Sudakov form factor Sa{Qf,Q^) denotes the 
probability for evolving from Qf to without splitting. 

We determine the virtuality of the parent parton 
according to the probability density i7a(£'^, Q^)- Consis- 
tent with the probability distributions Il^a^bc, we select 
the type of parton splitting, which occurred for parton a. 
The momentum fraction z of the splitting is then chosen 
within the kinematically allowed range 
z g [z^{Ql, E), z+i^Q"^, E)] for the decay of a parton of 
virtuality Q^, requiring that the virtualities of both par- 
tons are larger than the hadronisation scale parameter 
(3o/2- Subsequently, the virtual masses of the two daugh- 
ter partons are determined with the help of equation ([3]), 
subject to three constraints: The virtualities Qb, Qc of the 
daughters are required to be smaller than their energy z E 
or {l — z)E, respectively, and they must be larger than the 
hadronisation scale Qo/2. In addition, the virtual masses 
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of the daughters satisfy the constraint Q^+Qc < Qa- The 
branching a ^ 6 + c is finally accepted if the momentum 
fraction z chosen initially lies within the kinematically al- 
lowed range for daughters of 'mass' Qh and Qc, respec- 
tively. Otherwise, new values for Qb and Qc are chosen, 
in simulations with angular ordering it may also be nec- 
essary to reject the z value and try with a new one. If 
the branching is accepted, the full four-momentum is re- 
constructed for both daughters, assuming an azimuthally 
isotropic decay of the parent. The procedure is then iter- 
ated for the daughter partons, based on their assigned vir- 
tualitics. Parton evolution is stopped for daughter partons 
which do not split above the scale Qq and these partons 
arc declared to be on mass shell. Angular ordering is en- 
forced by allowing only splittings with decreasing emission 
angle. The coupling as is running at one loop. This shower 
is essentially the 'global' 'constrained' evolution which is 
one of the alternatives of the Pythia event generator. 

2.1.2 Hadronisation mechanism 

In order to have a hadronisation mechanism which is flexi- 
ble enough to also be useable for the complex parton state 
in heavy ion collisions, we have developed a modified ap- 
plication of the string hadronisation approach. The hadro- 
nisation prescription assumes maximal colour correlation 
between partons close in momentum space. The code iden- 
tifies first the parton with the highest energy in the event 
(alternatively, the highest p± may be used in hadron col- 
lisions.) If this parton is a gluon, it is split into a coUinear 
quark - antiquark pair with the energy sharing given by 
Pg^qq. The more energetic of the two is then the cndpoint 
of the first string. In the next step it is connected to the 
closest parton in momentum space (with the only excep- 
tion that a quark-antiquark pair from a single gluon split- 
ting is not allowed to recombine into a colour singlet). In 
case the closest parton is a gluon the procedure is iterated 
until either another (anti) quark is connected or there is 
no parton in the same hemisphere left. In the latter case 
a suitable cndpoint is generated by adding the required 
quark or antiquark with momentum in the beam direc- 
tion to the event. The whole procedure is repeated until 
all partons are connected in strings. The strings are then 
hadronised using the Lund string fragmentation [41] rou- 
tine of Pythia 6.4, with default values of hadronisation 
parameters. We checked that the strings are sufficiently 
massive for this string fragmentation routine to apply. 

This approach is inspired by the fragmentation of jets 
in hadronic collisions where the additional endpoints can 
be thought of as being part of the proton remnants, to 
which the jet is connected by colour fiow. The hadrons 
associated with this additional parton endpoint tend to 
go along the beam direction so that they are well sepa- 
rated from the jet. As long as the jet structure is analysed 
in a restricted rapidity range of approximately \Ari\ < 1 
around the rapidity of the parent parton, the resulting 
dependence of the model on this endpoint is negligible. 
There is the possibility to set a maximum invariant mass 
of neighboring partons in a string, that could be used to 
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Figure 1. The thrust, thrust major and thrust minor (Tr = 
(T, Tmaj , Tmin)) distributions for = 200 GcV e+e^ -^qq^ 
X collisions. Data of the ALEPH Collaboration [42] are com- 
pared to simulations of JEWEL: i) parton level after parton 
shower evolved down to Qo = 1 GcV, ii) hadron level after par- 
ton shower evolution followed by hadronisation {Qq = 1 GeV). 

tune the routine to data. However, we did not attempt to 

fine-tune this simplified hadronisation routine. 

This hadronisation mechanism has the advantage of 
being very flexible. It can be applied to any number of 
jets in an event and, more importantly, it can be applied 
also to jets in a nuclear environment. In the latter case, it 
may be desirable to hadronise only some high energy part 
of the event. This is also possible, since this hadronisa- 
tion routine can also handle systems that are not globally 
colour neutral. 

For e'^e" collisions, beam remnants do not exist and 
the two hemispheres are colour connected and do not 
hadronise separately. We therefore implemented a vari- 
ation of the string finding algorithm, which takes into ac- 
count these additional constraints on colour flow. For the 
observables discussed in this paper, we did not identify 
any significant difference between the two hadronisation 
mechanisms (MC results not shown). To simplify our pre- 
sentation, we therefore decided to simulate all observables 
with the same hadronisation routine described above, ir- 
respective of whether they are for e+e" or hadronic colli- 
sions. 

We checked that the hadronic observables are essen- 
tially independent of the cut-off scale Qq, at which the 
perturbative evolution is interfaced with string hadroni- 
sation. The default choice is Qo = 1 GeV. 



2.2 Comparison to data 

In this subsection, we compare the final hadronic state 
from Jewel to data on jets measured in y/s = 200 GeV 
e+e~ colHsions at LEP by the ALEPH collaboration [42]. 
To select the process e+e" qq ^ X at y/s = 200 GeV, 
these ALEPH data were taken with a veto on initial state 
radiation. They have been compared already [42] very 
favourably to standard event generators such as Pythia 
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Figure 2. The jet rates as a function of jet resolution scale 2/cut in — 200 GeV e+e~ ^ qq X collisions. Left hand side: 
Simulation of Jewel with and without hadronisation for evolution down to Qo = 1 GeV. Right hand side: Data of the Aleph 
collaboration [42] compared to simulations of Jewel with hadronisation. 



and Herwig. The purpose of the present data comparison 
is to validate the Jewel parton shower in the absence of 
medium effects against a set of benchmark data on jets 
measured at LEP, before studying the extension of the 
code to medium effects. 

To account for a 2-jet event. Jewel evolves each par- 
ent parton separately. For e~^e~ data, the code determines 
the initial virtuality of one of the two parent quarks ac- 
cording to the probability distribution Sg,{E^ ,Q^), given 
in equation ([3]). Then, there are two options: Either, the 
virtuality of the second parent quark is fixed to be the 
same, thus conserving energy and momentum at the hard 
vertex, or we choose it independently according to the 
same distribution ([3]). Since the differences observed be- 
tween these schemes turned out to be negligible, we present 
here only results for the latter case. 

We focus on three classes of tests. First, we test Jewel 
on the level of the overall energy flow, which is charac- 
terised by event shapes. Then we study an observable 
which is particularly sensitive to the parton fragmenta- 
tion pattern. Finally, we turn to a comparison of inclusive 
charged hadron distributions inside a jet. Hadronisation 
will be seen to play a more important role for the latter 
observables. 

In figure [H we compare data to simulations of three 
event shape observables, namely thrust T, thrust major 
Tinaj and thrust minor Tmin- For these quantities, one sums 
over the three-momenta pi of all final state particles. Ac- 
cording to the definition of thrust. 



on this direction n. 



T — maxn 



Ei \Pi ■ nrl 



(5) 



a 2-jet event is pencil-like if T = 1, that is if all particles 
are aligned parallel or antiparallel to a thrust axis rix- The 
event is spherical if T = 1/2. Once the thrust axis tit is 
known, one can determine the direction n orthogonal to 
riT, along which the momentum flow is maximal. Thrust 
major is defined as the projection of all particle momenta 



maj 



max„^..n.=o 



n 



(6) 



Thrust minor sums up the components of the final 
particle momenta Pi, which are orthogonal to the plane 
defined by n and ut, 



T,^ \Pi 



(7) 



We note that T, Tmaj and Tmin are perturbatively cal- 
culable, infrared-safe quantities and therefore particularly 
suitable for testing our parton shower implementation. 

As seen in figure [1] the final state parton shower pro- 
vides a reasonable description of these jet event shapes 
over most of the measured range. We recall that Jewel 
does not contain a matching of the parton cascade to ex- 
act matrix elements, which could improve the QCD mod- 
elling of large angle radiation. This may be the reason 
why the simulation gives fewer events with large Tmaj and 
Tmin- Also, we observe small deviations between data and 
simulation for very small values of Tmaj , T^in < 0.05 at 
the parton level. This region is known [43, 44] to be very 
sensitive to non-perturbative effects and hence cannot be 
reliably described by a leading order parton cascade as 
illustrated by the dashed curve. Consistent with this gen- 
eral statement, we find that for Tinaj , Tmin < 0.05, the 
inclusion of the hadronisation prescription leads to some 
improvement of the data comparison. 

We have compared Jewel to measurements of a wider 
class of event shape observables, in particular to oblate- 
ness, sphericity, planarity, aplanarity and total jet broad- 
ening. The comparison to these event shape observables 
are of similar or better quality than the comparisons shown 
in figure [1] With these studies, we have established that 
Jewel accounts for global features of jet energy flow with 
an accuracy which is sufficient to characterise (sufficiently 
large) medium effects on top of it. 
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Figure 3. The inclusive distribution dA'ch/d^, £, — 
In [-Ejot/Phadron] of charged hadrons in e'^e'^ qq X events 
at = 200 GeV. Data of the ALEPH Collaboration [42] are 
compared to simulations of Jewel: i) parton level after par- 
ton shower evolved down to Qo = 1 GeV, ii) hadron level after 
parton shower evolution to Qo = 1 GeV followed by hadroni- 
sation. 

In comparison to jet event shapes, there are measure- 
ments which are more sensitive to the discrete and stochas- 
tic nature of partonic processes underlying the QCD jet 
fragmentation. One such measurement is the so-called n- 
jet-fraction. Its definition is based on the Durham cluster- 
ing algorithm [45]. For each pair of final state particles, 
one defines a distance 

= 2min(i?f, E^){1 ~ cos%)/i?L ■ (8) 

The pair of particles with smallest yij is then replaced 
by a pseudo-particle, whose energy and momentum are 
the sums of its daughters. The clustering procedure is re- 
peated until all Uij exceed a given threshold j/cut- The 
number of clusters separated by a distance larger than 
Ucut is defined to be the number n of jets. Thus, as one 
decreases the resolution scale ?/cut, one becomes sensitive 
to finer and finer details of the discrete QCD radiation 
pattern. 

Figure [2] shows simulation results for the n-jet frac- 
tion. We find that for log]^Q(ycut) ~ — 3, the jet resolution 
scale is sufficiently coarse such that hadronisation plays a 
negligible role in the n-jet-fraction. It is only for smaller 
scales, that the hadronic late stage of QCD fragmenta- 
tion affects the number of jets identified by the Durham 
clustering algorithm. The 3-jet fraction is somewhat too 
small at large log]^Q(2/cut) while the 2-jet fraction is cor- 
respondingly too large, which is again due to the missing 
3-jet matrix element. We find that the n-jet fractions vary 
mildly with the size of the strong coupling on the partonic 
as well as on the hadronic level. Our choice of Aqqd leads 
to a reasonably good description of the data. 

In contrast to the measurements discussed so far, the 
modelling of single inclusive intra-jet hadron distributions 
and multi-hadron correlations requires detailed knowledge 
about the hadronisation mechanism. This is seen for in- 
stance in figure [3l where we compare results of our simu- 



lation to data of the inclusive distribution diVch/d^, { £, — 
In [Ejct/phadron]) of charged hadrons in e~^e~ ^ qq ^ X 
events at ^/s — 200 GeV. Irrespective of the scale Qo down 
to which the parton cascade is evolved prior to hadroni- 
sation, there is a marked difference between the hadronic 
and the partonic distribution. In particular, the partonic 
distribution is very 'hard', i.e. the yield of high momen- 
tum (small ^ < 1) partons exceeds the observed hadronic 
yield by far. Since high momentum partons are correlated 
in colour to softer partons, it is likely that any colour neu- 
tralisation mechanism, i.e. hadronisation, which takes into 
account colour flow, will fill the momentum space between 
these partonic partners. As a consequence, hadronisation 
is expected to soften the distribution for ^ < 1 consid- 
erable. This is seen for the string hadronisation mecha- 
nism, which accounts for the inclusive hadron distribution 
diVch/d^ (figure [3]) to a level better than a few %. 

We also made simulations with the independent string 
hadronisation mechanism [46-48], which does not invoke 
any colour correlation between partons and instead hadro- 
nises each parton in the final state separately. This mech- 
anism was found to depend strongly on Qq and tends to 
produce too many soft hadrons. 

We finally note that in the presence of a high-multi- 
plicity environment, novel hadronisation mechanisms may 
play a role. For instance, hadron formation may occur 
via recombination of partons [49-53] . Studying the hadro- 
chemistry of jets [35, 54] is likely to help characterising 
such novel hadronisation mechanisms, but a study of this 
is beyond the scope of this paper. 

Within these uncertainties, which are mainly related to 
the modelling of hadronisation, we have established that 
Jewel provides a reliable baseline for the characterisation 
of jet quenching phenomena. 

3 'Jet Quenching' in the Monte Carlo Model 

A Monte Carlo model of jet quenching requires a picture 
of how a final state parton shower is embedded in the 
space-time geometry of a nuclear collision, and how it in- 
teracts with the medium. The question where the par- 
ton splittings occur in space is likely to affect the result- 
ing medium modifications strongly, since it determines 
which partonic components interact for how long with 
the medium. However, very little is known about the spa- 
tiotemporal evolution of parton showers in the vacuum or 
in the medium. Inverting the standard logic, one may even 
view the medium as a probe of the spatiotemporal evo- 
lution of the jet, since the medium interferes with parton 
fragmentation on length scales comparable to time-dilated 
hadronisation times. This highlights that understanding 
the medium-modification of the jet relies on understand- 
ing how the medium interacts with it, and understanding 
properties of the medium relies on how they are reflected 
in the medium-modification of jets. Progress on this mu- 
tual dependence requires a dynamical description of both 
the jet and the medium. 

On general grounds, one expects that a parton of vir- 
tuality Q branches on a time scale 1/|Q| in its rest frame, 
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which translates to a Lorentz time-dilated length scale in 
the rest frame of the nuclear matter. Our model of the 
spatiotemporal distribution of parton branchings within a 
nuclear environment will be based on this relation between 
time scales and virtuality. While different parton showers 
use slightly different evolution variables, we have chosen 
virtuality (consistent with the Pythia 6.4 mass-ordered 
shower), since it makes the contact with a spatiotemporal 
picture relatively easy. 



3.1 Modelling the interactions with the medium 

We regard the medium as a collection of partons acting as 
scattering centres. For the case of elastic interactions be- 
tween the jet and the medium, each scattering centre dis- 
plays to the partonic projectiles an elastic 2-^2 scattering 
cross section da/dt. The form of this cross section, as well 
as the momentum distribution, the mass and the density 
of scattering centres can be specified freely in Jewel. For 
the studies presented in this paper, we focus on a specific 
model, in which the density and momentum distributions 
of scattering centres are that of a gas of massive quarks 
and gluons of temperature T. The masses of the scattering 
centres are fixed to mgcatt = fJ-D{T)/V2, where fJ.D{T) is 
the thermal Debyc mass. It is approximated by /id — 3T 
in the simulation. This model is minimal in the sense that 
the medium is characterised fully by a single parameter, 
the temperature T. 

To specify the spatiotemporal structure of the parton 
shower, we start from the estimate that the parton shower 
evolves down to components of virtuality Q/ on a time 
scale ~ For a parton of energy E and mass Q/, this 

will be time-dilated in the rest frame of the medium to 
~ E/Q'j. If the parton originated from the branching of 
some parton of virtuality Qi, then the parton of virtuality 
Q f existed for a duration of approximately 



E 

0} 



E 
01 



(9) 



For the parent parton, which initialised the parton shower, 
the lifetime is r = . In the case of a medium of constant 
density n, the probability that no scattering occurs during 
this time is given by 



S'noscattlT) = CXp [-CTdas " t/3] 



(10) 



Here, a^us is the total elastic scattering cross section and (3 
the projectile velocity. For the evaluation of the total scat- 
tering cross section in equation (jlOp the scattering centres 
are assumed to be at rest, for the simulation of the scat- 
tering process they are assigned a thermal momentum. If 
the density of scattering centres varies along the trajec- 
tory of the parton, one would replace this expression by 
an integral over the parton trajectory between initial and 

final times, exp — JJ^ crcias(0 "-(Od^ ■ Here, we have in- 
cluded the possibility that the elastic cross section changes 
during evolution, for instance, since its infrared regulation 



depends on properties of the medium, or as has a temper- 
ature dependence. 

For the differential elastic scattering cross sections, we 
choose 



d\t\ 



'-Cr- 



(11) 



regularised 



This is the leading i-channel exchange term for quark- 
quark {Cr = 4/9), quark-gluon {Cr = 1) and gluon-gluon 
{Cr = 9/4) scattering, respectively. The regularisation of 
equation (llip and the running of the coupling will be spec- 
ified below. Instead of equation (jlip . one could include the 
full LO 2 — > 2 parton cross sections [55]. These contain, 
for instance, also terms singular in the Mandelstam vari- 
able u. For a scattering in the medium, this corresponds 
to processes in which almost all the projectile energy is 
transferred to the scattering partner. Such processes may 
be viewed not as an energy degradation of the leading par- 
ton but as an exchange of the roles of projectile parton and 
target parton. For the sake of simplicity, the present study 
does not aim at including and characterising such more 
detailed features of elastic interactions between projectile 
and target. 

In the Monte Carlo model studied here, we know the 
virtuality of the partons, which emerge from parton split- 
tings or from elastic scatterings: After a splitting, the vir- 
tuality is determined via equation ([3]), after a scattering, 
it remains unchanged. This prescription is consistent with 
the dominance of small angle i-channel scatterings, which 
do not open phase space that could be used to reduce 
the fast parton's virtuality. Given this virtuality, we de- 
termine the lifetime r of the parton according to equa- 
tion With the probability 5'noscatt('''), this parton will 
not scatter, and the code simulates the next splitting as in 
the vacuum case. With the probability [1 — S'noscatt(T)], 
the parton scatters at some time r' < r. In this case, 
the parton exchanges momentum with a scattering centre 
according to the differential cross section equation pT|). 
Then, the code continues to propagate the scattered pro- 
jectile parton, by checking whether there will be further 
parton scatterings before splitting. In this procedure, an- 
gular ordering is enforced for consecutive splittings, but 
it is reset at each scattering, consistent with the assump- 
tion that scattering destroys the interference pattern of 
parton radiation. On-shell partons are allowed to scatter 
until they leave the medium. It is thus assumed that par- 
tons cannot hadronise inside the medium. The code has an 
option to propagate also the recoiling target partons with 
subsequent splittings and scatterings. To simplify the pre- 
sentation and discussion, we do not explore this option in 
the present work. 

The above prescription for embedding a parton shower 
in a nuclear environment includes inevitably model-depen- 
dent assumptions. In this exploratory work, we do not em- 
bark on their systematic study, but we would like to recall 
some of the major sources of such model dependencies: 

First, as mentioned already, the spatiotemporal struc- 
ture of the parton shower is based on assumptions, which 
are difficult to constrain. For instance, even if the lifetime 
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r of a virtual state is given parametrically by ([9]), one may 
distribute this lifetime with an exponential decay law with 
half time ([9]), or one may try to improve the lifetime es- 
timate e.g. by accounting more accurately for the energy 
of the daughter partons or by replacing it with other as- 
sumptions. 

Second, the choice of the cross section ([TT|) invokes the 
assumption that the interaction between the jet and the 
medium can be treated perturbatively. This assumption is 
problematic, since data at RHIC provide some evidence of 
anomalously strong coupling. In addition, there are quan- 
titative uncertainties arising from the regularisation of the 
cross section (fTTj) . 

Third, the model introduced here does not yet include 
a mechanism of radiative energy loss. Since this may be 
the main source of energy degradation, we have included 
an option to enhance the vacuum splitting functions by a 
factor (1 + /mod), 



Pa 



(l + /n 



(12) 
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as long as the splitting occurs in the medium. This pre- 
scription has been argued to mimic characteristics of ra- 
diative energy loss [34] . In a probabilistically iterated par- 
ton shower, a more satisfactory treatment of radiative en- 
ergy loss may be to include partonic inelastic 2 — > 3 scat- 
tering cross sections on the same footing as the elastic 
2^2 cross sections included here. We did not imple- 
ment this option in the present study, since its discussion 
involves significant additional conceptual and technical is- 
sues, which - in our view - deserve a separate study. 

Fourth, there is the question of how to interface a 
medium-modified parton shower with a hadronisation model 
In the present study, the final partonic state is interfaced 
with the hadronisation model discussed in section [2] after 
all interactions with the medium. So, the present version 
of Jewel is based on the assumption that the hadroni- 
sation of a jet fragment is unaffected by medium effects, 
but this assumption could be modified in future studies to 
explore for instance changes in jet hadrochemistry [35]. 

The present paper does not aim at fully studying or 
resolving any of these issues. But we expect that Jewel 
provides a tool for their systematic exploration in subse- 
quent works. 



Figure 4. The average parton energy loss d_E/da; of a quark 
of energy E, undergoing multiple elastic collisions over a path 
length L — 1 fm in a thermal medium of temperature T. 
Elastic collisions are described by the infra-red regulated par- 
tonic cross sections of equation (I13|l (case I) and equation (|14p 
(case II). 



scattering centres. Here, we only consider the coUisional 
energy loss of massless partons and we leave the case of 
heavy flavours to a dedicated study. 

As discussed in section l3Tl we define the density, mass 
and momentum distribution of scattering centres in terms 
of a thermal distribution. The same temperature is the 
model parameter entering the regularisation of elastic scat- 
tering cross sections. To make contact with the previ- 
ous studies, we haven chosen two different regularisation 
schemes for the elastic scattering cross section, namely 
case I 



|t.„ax| 



clas 



d\tr-^^Mp3Ac,-i±ii^, (13) 



which is the default, and case II 



^clas 



d\t 



naU\t\) 



Cr 



+ (^-1^1)' 



(14) 



3.2 A baseline of medium-efrects: collisional energy 
loss without parton branching 

The case of a single, high energy parton losing energy via 
multiple elastic interactions in a spatially extended tar- 
get has been studied extensively in the literature [8-16]. 
This problem can be studied in Jewel by switching off 
the option of parton splitting, and specifying the (model- 
dependent) density of scattering centres and elastic scat- 
tering cross sections. The MC then simulates the propa- 
gation of an on-shell parton (massless or massive) which 
undergoes scattering in a medium. The scattering proba- 
bility is in this case determined only by the path length 
inside the medium, and the density and cross section of 



The kinematic boundary is [tmaxj = 3^/2ET. The total 
cross sections (fT3|) and (|14|) differ only by their regularisa- 
tion. In figure 01 we have calculated the resulting average 
energy loss for an in-medium path length of L = 1 fm in 
a medium of temperature T. In general, the average col- 
lisional energy loss increases with increasing temperature 
and projectile energy. Quantitatively, the models I and II 
show differences of approximately a factor 2 in dE/dx for 
a 10 GeV parton. The differences decrease slowly with in- 
creasing projectile energy. Cross section II leads to a larger 
energy loss, as may be expected since there is minimum 
momentum transfer. 

We note that in calculations based on finite tempera- 
ture field theory, the medium specifies a preferred Lorentz 
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Figure 5. Distribution of energy loss AE for different parameter choices (top panels) and the probability for an energy loss 
smaller than AE after passage of a medium of length L (bottom panels) . 



frame and the effect of collisional energy loss is in general 
not fully described by an elastic scattering cross section 
of the form da/dt. This complicates analytical compar- 
isons of our model with recent studies of collisional energy 
loss [8-16]. On the other hand, to the best of our knowl- 
edge, no characteristic medium-induced deviations from 
the form dcr/dt have been explored so far, and a compar- 
ison of our calculations for dE/dx on the numerical level 
seem meaningful. We find that the temperature depen- 
dence of d£'/da; shown in figure [H is consistent with the 
dependences reported previously. The differences between 
case I and case II are representative of the typical factor 
2 uncertainties between different model studies, though 
some recent studies lead to slightly larger values of dE/dx 
than those shown in figure [4j see e.g. [14]. Here, we do 
not enter the current discussion of whether an improved 
understanding of the scale dependence of elastic interac- 
tions may allow to narrow these uncertainties. Rather, we 
take equation p3p as an example of an elastic cross sec- 
tion, whose strength and t-dependence may be changed 
in future simulations of Jewel. We finally note that in 
a Monte Carlo simulation, the average energy loss does 
not grow linearly with the in-medium path length. This 
is because energy and momentum are conserved exactly 
so that subsequent collisions occur at decreasing centre of 
mass energies, which introduces a mild non-linearity. 

For single inclusive spectra, which fall steeply in trans- 
verse momentum, the dependence of parton energy loss is 



not determined by the average energy loss but by the least 
energy loss which a significant fraction of the projectiles 
suffer. For this reason, it has become the standard in stud- 
ies of collisional energy loss to study dN/ dAE of partons 
suffering an energy loss AE. Our result for this distri- 
bution is shown in figure [5] and is seen to peak sharply 
for negligible energy loss. This refiects the fact that elas- 
tic interactions are dominated by small angle scattering, 
for which the longitudinal component of the momentum 
transfer between projectile and target becomes negligible 
in the high energy limit. 

We note that the temperature dependence of colli- 
sional parton energy loss models shows several generic 
features. Qualitatively, these may be understood by ob- 
serving that the strength of the average collisional energy 
loss is essentially governed by the product of the density 
of scattering centers times the cross section, which each of 
these scattering centers displays. The density rises propor- 
tional to with increasing temperature. The cross sec- 
tion, however, is dominated by the lowest typical momen- 
tum transfer, which is governed by the Debye screening 
mass, so parametrically a cx ^/m^ cx For increas- 

ing temperature, the strong increase in density is largely 
balanced by a reduction of the effective cross section per 
scattering center, and the net effect is an average energy 
loss AE which shows a mild, approximately linear increase 
with T, as seen in Fig. [H For the same reason, increas- 
ing the Debye screening mass by increasing T, one re- 
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Figure 6. Transverse momentum relative to beam axis and angle with respect to jet axis of recoiling scattering centres as 
compared to the undisturbed medium for different temperatures, with and without splitting of the projectile (cross section I, 
Ejet = 100 GeV, L = Ifm). The jet is at mid-rapidity (?7 = 0). Hadronisation is not included but may affect these distributions 
significantly (see text for further discussion) . 



duces the relative yield of scatterings at small t strongly, 
while the rate of particles propagating without interac- 
tion through a finite size medium is affected only mildely. 
As a consequence, the peak in the distribution dN/dAE 
at AE = depends on details of the infrared regulariza- 
tion of the elastic cross section, but tends to become more 
pronounced with increasing temperature. The mean free 
path changes only mildly with temperature but is rela- 
tively large (1.7 fm for a 100 GeV quark at T = 500 MeV). 

Another characteristic property of t-channel dominated 
elastic differential cross sections of the form equation (jlip 
is their weak y^-dependence in the high energy limit. This 
leads to a remarkably weak dependence of dN/dAE on 
projectile energy E, as seen in figure O We note, how- 
ever, that the dN/dAE distribution extends up to the 
limit AE = E. As a. consequence, the average energy loss 

= lo <ii^E) ^ AE / J^^ d{AE) ^ differs sig- 
nificantly for E = 10 GeV and E = 100 GeV, since for 
larger jet energy E, the integrand AE has contribu- 
tions for larger values of AE. As seen in figure SI there is 
a marked £'-dependence of the average energy loss, while 
there is a very weak £'-dependence of for small val- 
ues of AE. This illustrates the well-known fact (for more 
details, see e.g. Ref. [26]) that the value of the average en- 
ergy loss is dominated by rare events with large AE, and 



thus —dE/dx does not characterise adequately the 'typi- 
cal' energy loss suffered by most events entering -^g- 

Upon proper normalisation, the distribution dN/ dAE 
can be turned into a probability distribution for a high 
energy quark to lose less than a specific amount of energy 
AE. This quantity is plotted in the lower panel of figure [5] 
One sees that for an in-medium path- length of L = 1 fm, 
almost 65 % of all projectile partons do not lose any en- 
ergy in our models, even if the temperature is taken to be 
500 MeV. For 3fm in-medium path length, there are still 
typically 25 % of all projectiles which emerge unscathed. 
This indicates that elastic interactions alone are unlikely 
to account for a large fraction of the suppression of single 
inclusive hadron spectra. 



3.3 Characterising the Recoiling IVIedium 

The energy lost by a jet is redistributed amongst the tar- 
get components. If one could characterise the amount of 
energy recoiling in the medium, this would be an unam- 
biguous determination of coUisional effects. (We note that 
also 2^3 and 2 — s- n processes can contribute to such 
recoil effects.) More generally, characterising the recoil of 
the medium may provide a means to disentangle differ- 
ent mechanisms of parton energy loss. In this section, we 
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Figure 7. Correlation of angle and momentum of recoil- 
ing scattering centres in the Monte Carlo simulation (T = 
500 MeV, without splitting, L = 1 fm, £ = 100 GeV) and ana- 
lytical estimate (equation p8|l for mean Bin) 



characterise the momentum distribution of target partons, 
recoiling against a projectile parton which undergoes elas- 
tic 2 —> 2 processes in a thermal medium. 

In figure [6] the momentum and angular distribution 
of recoiling scattering centres is compared to the undis- 
turbed thermal medium. We investigate the cases, that 
i) an on-shell projectile is propagated without splitting 
(upper panel of figure [5]) or that ii) the splitting of the 
projectile is included (lower panel of figure [5]) . The results 
are very similar for both cases, indicating that the recoil 
is only weakly sensitive to the projectile energy. 

The scattering centres in the medium show initially 
an exponentially falling p_L-spectrum and an isotropic mo- 
mentum distribution, characteristic for a thermal distribu- 
tion. After interaction with the projectile, thepx-spectrum 
of the recoil follows a power law at intermediate and high 
p^. The high-p^ tail is more pronounced at higher tem- 
peratures. Using the kinematics of 2 — > 2 scattering, one 
estimates 



l^l+2i^fn 

2-EiTi 



da „ d(T , , 

Mzr'-^'-'Wv '''' 



where E^^ and -Bout are the energies of the incoming and 
outgoing scattering centre, respectively, and the average 
over the angle between the initial partons was taken. We 
note that despite this power law, the yield of recoil par- 
ticles lies mainly at relatively soft momentum, p±^ < 1 — 
2 GeV. 



The recoil moves predominantly in the jet direction 
with a maximum at A(j) ~ 0.8 nearly independent of the 
temperature. The results shown here are for cross sec- 
tion I; with option II the result is qualitatively the same. 
There is a moderate change in the shape of the distribu- 
tion, but the position of the maximum remains unchanged. 

Most of the features of the angular correlation can be 
understood from a simplified analytical model where the 
scattering centre is assumed to be at rest and the cross 
section is approximated by dCT/d|t| (x (very similar 

results are obtained with da/d|i| oc {\t\ + /i^)"^). In this 
approximation cos(Z\(/<) can only be positive. Given the 
shape of the distribution the position of the maximum 
can be approximated by the expectation value 



(cos(Z\0)) ~ (\/3- 1 



™scatt 



2^3-3 



(16) 



which is practically constant for large projectile energies 
Ep. It has, in particular, only a weak temperature depen- 
dence through the mass mscatt of the scattering centre. 
However, the position of the maximum is determined by 
the regularisation of the cross section. We therefore inves- 
tigate the dependence on the infrared regulator by treat- 
ing it as an independent variable. After expansion around 
/.tp = 2mg(.a^tt one gets for the shift of the maximum posi- 
tion 



Zi(cos(zi0)) 



%/3- 1 



1 



2Tr)2 
^'"-scatt 



(17) 



This finding is in agreement with the results from the full 
Monte Carlo simulation. 

There is a strong correlation between the scattering 
angle and the momentum transfer, such that the most en- 
ergetic partons are closer to the jet axis. If one allows the 
incoming scattering centre to carry momentum but aver- 
ages over its direction, then the analytical model discussed 
above leads to an angular dependence of the outgoing scat- 
tering centre with respect to the incoming projectile, 



cos{A(j)) 



Eont {Ep + Ein) 



EpEu 



El 



(18) 



In figure[7]this result is compared to the Monte Carlo sim- 
ulation. The overall behaviour is similar, but the Monte 
Carlo result is distributed wider due to the spread in en- 
ergy and angle of the incoming scattering centres. 

At face value, figure [5] indicates that a jet can be ac- 
companied by additional multiplicity which has its maxi- 
mum separated from the jet axis by a characteristic finite 
angle Acj), and whose yield is expected to die out quickly 
with increasing transverse momentum. Similarly, figure [7] 
indicates that scatterings with large energy transfer lead 
to a relatively large projectile energy loss of the leading 
parton, but that the recoil parton is very close to the jet 
axis so that the energy may stay within the jet cone. Typ- 
ical momentum transfers, on the other hand, tend to scat- 
ter the recoil to a characteristic and relatively large angle. 
We note, however, that whether the partonic distributions 
of figure inland [7] will or will not change substantially upon 
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hadronisation may depend on details of the hadronisa- 
tion mechanism. In particular, if the partons entering the 
peak of the cos(Z\(/))- distribution are connected by colour 
strings with the leading jet fragments at cos(Z\0) ^ 1, 
then the hadronic distribution will peak much closer to 
cos(Z\0) = 1 than the partonic one shown in figure [B] On 
the other hand, multiple gluon exchanges with the target 
break the colour flow between hard projectile fragments 
and target recoils. If such a mechanism is invoked, then 
the hadronic distribution is likely to follow the partonic 
one shown in figure [51 Here, we limit our discussion to 
the kinematic constraints of the underlying partonic dis- 
tributions in figure [6] and [3 but we do not explore the 
model dependence associated with the hadronisation of 
these distributions. 



4 Medium-modifications of jet measurements 

In this section, we study to what extent the jet measure- 
ments discussed in section [2] are sensitive to the medium 
effects introduced in section [3l A medium does not only 
affect jet fragmentation. It also results in a high 'back- 
ground' multiplicity of soft particles. This complicates the 
characterisation of jets and their medium modifications. 
To establish which jet measurements are experimentally 
feasible despite this background, one would ideally like to 
embed the parton shower simulated by Jewel into full 
event simulation of a heavy ion collision. This could be 
done but it lies beyond the scope of the present paper. In 
this section, we compare analyses of the full simulated jets 
with analyses, for which we only include hadrons above a 
background cut of £^cut = 2 GeV. This comparison may 
indicate the extent to which different jet medium modifi- 
cations remain visible above the hadronic background of 
the heavy ion collision. 

For e'^e~-collisions, thrust, thrust major and thrust 
minor are of particular interest for testing QCD radiation 
physics, since they are infrared safe and thus perturba- 
tively calculable quantities. Applying the definitions (O, 
([U and (O to a single jet or to the jet activity above 
background, the feature of infra-red safety is lost. But T, 
Tmaj and Tmin will still provide a characterisation of the 
global energy flow in a jet. Here, we ask to what extent 
these quantities may provide a useful characterisation of 
jet medium modifications. Figure [8] shows the result of 
simulations of the Jewel parton shower. If coUisional en- 
ergy loss is included, then the jet is expected to broaden. 
This is clearly seen in the broadening of T, Tmaj and Tmin 
in the upper left panel of figure \8\ However, the kinemat- 
ics of elastic 2^2 scatterings dictates that the more 
energetic projectiles are deflected by smaller angles. We 
checked that the collisional broadening observed in fig- 
ure [8] is due to recoil. If the recoiling scattering centres 
are removed from the final state and do not hadronise, the 
thrust distributions with medium are practically indistin- 
guishable from the vacuum results. The recoiling scatter- 
ing centres have mostly relatively low momenta so that 
the medium-induced broadening is much reduced and be- 
comes small if only hadrons of energy above Ecut = 2 GeV 



are included in the analysis (see lower left panel of fig- 
ure[8|). In contrast to collisional energy loss, single partonic 
components of medium-induced additional radiation have 
a higher probability to carry a significant energy fraction 
of the initial projectile energy. So, on general grounds, one 
expects that the medium-induced broadening of the dis- 
tributions in thrust, thrust major and thrust minor will 
persist even if soft hadrons of energy £'h < Ecut = 2 GeV 
are dropped from the analysis. This is clearly seen in fig- 
ure [H We note that the value of /,ned entering these sim- 
ulations is a priori a free model parameter. We estimated, 
however, that the choice /mod = 3 lies within the range of 
parameters consistent with a nuclear modification factor 
RphPh ~ 0.2 in central lead-lead coUisions. 

The medium-modification of the n-jet fraction may 
provide another tool for disentangling elastic from inelas- 
tic jet quenching mechanisms. The potential impact of 
elastic 2^2 processes on the n-jet fraction is limited to 
two numerically rather minor effects: First, elastic colli- 
sions are a source of angular broadening, which can in- 
crease the distance ytj between partons and this implies 
that additional jets can become visible on a coarser scale 
Ucut- However, the higher the energy of the partonic pro- 
jectile, the smaller is this broadening effect. For this rea- 
son, one expects the dependence of the n-jet fraction on 
collisional energy loss to be relatively small and mainly 
visible for relatively small values of ?/cut which are sen- 
sitive to small energies. Second, collisional mechanisms 
may kick target components towards higher transverse 
momenta. In this way, hard recoils may be counted as 
jets, thereby increasing the n-jet fraction. However, from 
the momentum space distribution of recoil partners stud- 
ied in figures [6] and [71 we know that this effect will 
remain limited to very fine jet resolution scales. This is 
confirmed by the results shown in figure [3 which shows 
that even after an in- medium path length of i = 5 fm in 
a medium of temperature T — 500 MeV, the n-jet frac- 
tion remains unaffected above log]^o(ycut) = — 1-5 and the 
mild medium-induced deviations seen at higher jet reso- 
lution scale can be attributed largely to soft hadrons with 
energy Ei, < E^ut = 2 GeV. 

In contrast to elastic interactions, medium-induced ra- 
diation results in a distribution of sub-leading jet frag- 
ments up to characteristically larger energies. This is clearly 
seen in the medium-modified n-jet fractions of figure [9l 
which shows a marked increase in the average number of 
jets even for a very coarse resolution scale log]^o(ycut) > 
— 1. After a background cut at E^ut — 2 GeV, this medium- 
modification remains clearly visible on the logarithmic 
scale logio(2/ cut } • 

We have checked that the n-jet fraction shows a very 
weak dependence on the total jet energy in the range 
75 < i?jct < 125 GeV. However, if this jet energy is not 
measured exactly, E'jot.truo = /corr -Ejet ,meas, then the n- 
jet fraction is shifted by a constant value —\ogiQ{fco„) 
in log2^Q(ycut)- For instance, a deviation of the jet en- 
ergy reconstructed from a calorimetric jet measurement 
from the true jet energy by 30% amounts to a shift of 
log]^Q(1.3^) « 0.23. This is comparable to the size of the 
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Figure 8. Thrust, thrust major and thrust minor (Tr — (T, Tmaj , Tmin)) for a single 100 GeV jet. The Jewel parton shower 
in the vacuum is compared to two scenarios including medium-induced parton energy loss. Left hand side: Collisional energy 
loss for a medium of T = 500 Me V and in-medium path length L = 5fm (the recoil is hadronised with the medium). Right 
hand side: Radiative energy loss for /med = 3 and L = 5fm. For the plots in the lower panel, only hadrons with energy above 
_Ecut ~ 2 GeV are included. 



medium modifications on the right hand side of figure [9l 
So, our studies indicate that the n-jet fraction may provide 
a powerful tool for disentangling radiative from collisional 
energy loss mechanisms, but a proper assessment of var- 
ious experimental uncertainties will be clearly needed to 
explore this tool. An analogous statement applies to the 
medium-modifications of jet shapes studied in figure [H 

In figure IIO) we plot the medium modification of the 
single inclusive distribution dN/ d^ in the presence of col- 
lisional and radiative medium effects. Similar to the case 
for the vacuum shower in figure [3l hadronisation is seen to 
lead to a significant softening of the distribution. If elas- 
tic interactions with the medium are included but only 
the cascade is hadronised, then the total jet multiplicity 
of projectile partons or hadrons depends only very weakly 
on the medium, since 2^2 processes do not increase 
the parton multiplicity. However, the total jet multiplic- 
ity may increase significantly if recoil partons are counted 
towards the jet. This is illustrated in figure fTOl which also 
shows the resulting dN/d^ distribution after addition of 
all recoil partons that go in the same hemisphere as the 
original parton. Thus the jet-modification of the medium 
and the medium-modification of the jet are two comple- 
mentary aspects of the same dynamic phenomenon. 



As stated above, the present study of radiative parton 
energy loss mechanisms is based on a simplified model 
which does not account for the recoil distribution of tar- 
get partons. However, figure [10] indicates that radiative 
mechanisms result in a larger increase of intra-jet multi- 
plicities. 



5 Conclusion and Outlook 

Monte Carlo models for hadronic and nuclear collisions 
are at the interface between the theory of QCD and ex- 
periment. Depending on the class of measurements and 
the status of theory, they provide a bridge between QCD 
and data which invokes a varying degree of model as- 
sumptions. In this paper, we have started to discuss the 
Monte Carlo model Jewel, developed for a class of mea- 
surements sensitive to jet quenching phenomena in ultra- 
relativistic nuclear collisions. In the absence of medium- 
effects, we have validated Jewel in section [2] against a 
set of benchmark jet measurements. In section [31 we have 
then introduced medium-modifications with a special fo- 
cus on the dynamic description of elastic interactions be- 
tween a fragmenting projectile and the medium. The main 
results, presented in section [4l show that collisional energy 
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Figure 9. n-jet fractions for a single 100 GeV quark jet after hadronisation in vacuum and with medium effects. Left hand side: 
Collisional energy loss for T = 500 MeV and L = 5fm. Recoil partons are either hadronised together with the cascade ('all') or 
they are not included in the hadronisation ('cascade'). Right hand side: Medium-induced radiation for fmcd ~ 3 and L = 5fm. 
In the top row, all hadrons are included, while in the bottom row, only hadrons with energy above _Ecut ~ 2 GeV are included. 



loss is important and demonstrate the basic physical jet 
quenching effects for given temperature and path length 
in the plasma. Moreover, we find that quantities like the 
n-jet fraction and jet event shapes may allow us to dis- 
tinguish between elastic and inelastic mechanisms under- 
lying medium- induced parton energy loss. As discussed 
in the introduction, the microscopic mechanism(s) under- 
lying jet quenching are not firmly established and their 
theoretical description is incomplete. In view of the com- 
plexity of this problem in heavy ion collisions, the Monte 
Carlo simulation method is particularly suitable for a de- 
tailed treatment. Amongst the many open problems in the 
theory and phenomenology of jet quenching, we identify 
in particular the following points which may be accessible 
by further developments of Jewel: 

1. Heavy flavour, multi-hadron correlations and other mea- 
surements. 

Jewel can be extended easily to include the in-medium 
propagation of heavy quarks, which has attracted sig- 
nificant interest recently [11-14,56]. Simulations of multi- 
hadron correlations can be extracted from the present 
version of Jewel, but are likely to depend strongly on 
the hadronisation model, as discussed in the context 
of figure m 



2. Realistic description of collision geometry and expan- 
sion. 

On a schematic level, Jewel can be interfaced with 
any simulation of the bulk evolution of heavy ion colli- 
sions, by supplementing Jewel with information about 
a realistic distribution of partonic trajectories inside 
the medium, and with information about a realistic 
time-dependent density of scattering centres along these 
trajectories. In particular, this includes full 2-dimen- 
sional ideal [57] or viscous [58,59] and full 3-dimensio- 
nal ideal [60] hydrodynamic simulations of the heavy 
ion collision, as well as other models for an expanding 
plasma [15]. 

3. Radiative parton energy loss. 

The version of Jewel presented here mimics radia- 
tive energy loss by enhancing the pcrturbativc split- 
ting functions. We regard the substitution of this ad 
hoc prescription by realistic 2 — > 3 inelastic parton 
scattering cross sections as a feasible and interesting 
future extension of Jewel. It would provide for a ra- 
diative energy loss mechanism which includes recoil 
effects, and it could define the incoherent limit of any 
radiative mechanism. 

4. The separation between weakly coupled and strongly 
coupled regimes and the problem of hadronisation. 
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On theoretical grounds, there is ample motivation to 
use weak coupling techniques for describing the split- 
ting of highly virtual partons, and their sufficiently 
hard interactions with a medium. However, for soft in- 
teractions with the medium, one has to invoke a non- 
perturbative mechanism for regularising the infra-red 
singularities, resulting in uncertainties discussed in sec- 
tion [3l Additional complications arise if one considers 
the propagation of soft projectile fragments whose mo- 
mentum is close to that of components of the heat 
bath. Propagation of these fragments within a pertur- 
bative parton cascade is expected to be unreliable since 
the medium is likely to be strongly coupled and the 
dynamics of the medium cannot be accounted for by 
partonic 2^2 and 2^3 processes only. Another 
open question is how to hadronise these soft jet com- 
ponents within a soft high multiplicity environment, 
where novel hadronisation mechanisms such as recom- 
bination may become relevant. These processes may 
become a test laboratory for understanding how a well- 
defined partonic projectile interacts and to what extent 
it thermalises within a finite size medium. From a prag- 
matic point of view, a parton shower may contribute 
to this issue by identifying the momentum scales at 
which a perturbative description breaks down. 

There are many other open issues for improving our 
understanding of jet quenching in an interplay between 
theory and experiment. One of the most central ones may 
be how to best characterise a jet within a high multi- 
plicity environment such that unambiguous information 
about its medium modifications can be obtained. Stan- 
dard calorimetric jet measurements are very difficult to 
apply to heavy ion collisions, though novel jet algorithms 
may be better suited [61]. On the other hand, the suppres- 
sion of single inclusive spectra provides unambiguous evi- 
dence that there are strong medium modifications. There 
is no fundamental reason why a similar unambiguous char- 



acterisation should not be possible for other aspects of the 
entire jet fragmentation pattern. From a pragmatic point 
of view, the question is then how to identify classes of jet 
measurements, which are sensitive to medium effects but 
which remain sufficiently insensitive to operational uncer- 
tainties in the jet definition. Jewel may contribute to 
this central issue on various levels. By superimposing sim- 
ulations of Jewel on top of the simulated background 
of heavy ion collisions (or developing Jewel to include 
this) one can test the sensitivity of different jet observ- 
ables. On the other hand, by simulating the redistribution 
of 'background' multiplicity due to the propagation of a 
jet. Jewel provides a mean to go beyond the simplifying 
assumption that the medium-modified jet is uncorrelated 
to the underlying background. These features of Jewel 
may contribute to establish to what extent an operational 
procedure of characterising jet medium modifications is 
suited to draw model-independent conclusions. 
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